function dXdt = inductive(t, X, param)
% 感应耦合
% 状态变量 X = [x1; y1; x2; y2; z]

    x1 = X(1);
    y1 = X(2);
    x2 = X(3);
    y2 = X(4);
    z  = X(5);

    b = param.b;
    r = param.r;
    c = param.c;
    n = param.n;
    lambda_L = param.lambda_L;
    
    % 混沌电流和外部声信号
    I_Chaos = param.iChaosFunc(t);
    u_PC = param.vPCFunc(t);

    % 根据式(6)定义微分方程
    dx1 = x1 - (1/3)*x1^3 - y1 + I_Chaos + z;
    dy1 = c * (x1 + b - r*y1);
    dx2 = x2*(1 - eta) - (1/3)*x2^3 - y2 + n*u_PC - z;
    dy2 = c * (x2 + b - r*y2);
    dz = lambda_L * (x2 - x1);

    dXdt = [dx1; dy1; dx2; dy2; dz];
end